A Droplet State in an Interacting Two-Dimensional Electron System 
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It is well known that the dielectric constant of two-dimensional (2D) electron system goes negative 
at low electron densities. A consequence of the negative dielectric constant could be the formation 
of the droplet state. The droplet state is a two-phase coexistence region of high density liquid 
and low density "gas". In this paper, we carry out energetic calculations to study the stability 
of the droplet ground state. The possible relevance of the droplet state to recently observed 2D 
metal-insulator transition is also discussed. 
PACS numbers: 71.30.+h, 73.40. Hm 
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The recent discovery of two-dimensional (2D) metal- 
insulator transition (MIT) by Kravchenko et al. M has 
challenged the scaling theory of localization j^,|3j in which 
a 2D MIT is forbidden. A noticeable character of the elec- 
tron system in these experiments is that r s , the parame- 
ter measures the strength of the Coulomb interaction, is 
fairly large. We suspect that the electron system may be 
unstable against phase separation at these large values 
of r s . We demonstrated in our previous paper || that 
this assumption alone is sufficient to provide a theoreti- 
cal description that is consistent with all the known ex- 
periments. For a two-dimensional (2D) electron system, 
there believed to be two phases: a high density Fermi gas 
phase and a low density insulating Wigner crystal phase. 
The dielectric constant of the liquid phase becomes neg- 
ative when r s ~ 2 |f|, which indicates that the liquid 
phase is unstable. At lower densities, the Wigner crystal 
phase appears around r s ~ 37 in the absence of disorder 
||. This critical value of r s appears to be reduced to 
around r s ~ 10 with disorders Q. In the intermediate 
values of r s , we believe that there is a liquid phase which 
we think is responsible for the observed MIT. 

In this paper, we propose that a droplet state of the 
electron system resulted from the phase separation of the 
electrons into this new liquid phase and a low density 
"gas" phase. Here we call the low density phase "gas" 
purely for the reason that its density is low. In fact, 
in the presence of impurities, the "gas" phase is disor- 
dered Wigner crystal. To investigate our proposal, we 
have studied the energetics of such a droplet state. We 
find that both electron-electron interaction and potential 
fluctuations are crucial for the formation of the droplet 
state. 

An obvious condition for the droplet state is that the 
electron gas is unstable. To investigate what possibil- 
ities of the instability leads to, we study a simple but 
physically motivated model. Let us consider electrons in 
the disc of radius b with positive background. Imagine 
that the electron system is shrunk to a new radius a < b 
while the positive background remains intact. Clearly 
the charging energy due to the separation of the elec- 



trons from the positive background increases the energy 
of the system. However there can also be energy gained 
(decreasing total energy): since for a uniform electron 
gas the ground state energy E g is at its minimum when 
r s ~ 2 j^], for r s > 2 the system gains energy by shrink- 
ing the area occupied by electrons. Furthermore, in the 
presence of disorder, electrons tend to occupy the valleys 
of the disorder landscape. Thus, a slowly varying disor- 
der potential is in favor for the formation of the droplet 
state. We calculate the energy changes when electron 
disc is shrunk from b to a to determine whether a spon- 
taneous shrinking can take place. 

For electron-electron or electron and positive back- 
ground interaction, we use screened Coulomb potential. 
For Si MOSFETs, the image charge in the metal sub- 
strate induces the screening and the interaction in the 
momentum space can be written as [0 
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where D is the thickness of the S12O insulating layer 
and e — \ (ei +£2), K = (ei - e 2 )/(ei + £2), with ei 
and e 2 being the dielectric constants of Si and S^O, 
respectively. For other systems, such as GaAs/AlGaAs, 
the screened interaction can be well represented by the 
following form: 
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and the corresponding moment space representation is, 

V ' e V*: 2 + A 2 

Both forms of the interactions define an interaction range 
£. For Si MOSFETs, £ = —D, and for the screened 

Coulomb potential, S, — j\- Outside the range the inter- 
action is strongly screened. 
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FIG. 1. The ground state energy of electron system for 
different screening strength: (a) for screened Coulomb inter- 
action, (b) for Si MOSFETs. 

To calculate the ground state energy of the electron gas 
with the screened Coulomb interaction, we use the varia- 
tional correlated-basis-function (CBF) method j9|. This 
method has been applied to the bare Coulomb interac- 
tion and is proved to provide rather accurate result for 
the ground state energy of 2D electron system ||. The 
accuracy of the ground state energy is found to within 
10% comparing with the best available quantum Monte- 
Carlo results H for densities down to r s — 20. In the 
CBF approach, there is a variational variable a for which 
the ground state energy E g has to be minimized 
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where lRy = a B = and y u = y x - y 2 - 

U c ("7^) i s ^ ne cohesive energy which depends on the 
special form of the interaction, 
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Figure [l] shows the calculated results for the ground 
state energy. The screening effect raises the ground state 
energy because the electron correlation is suppressed. 

The shrinking of the electron disc will cause the re- 
distribution of the charge, which will raise extra electro- 
static energy because the positive background is fixed. 
The charge distribution can be written as 



p{k) = 2N 
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where N is the total number of the electrons in the disk. 
The charging energy is 



1 



d 2 k 



E c = 7^ T 77^V{k)\p{k)\ 



2N J (2tt 



1 



na^r 



B s JO 



dx /x 

~ V Kb 



Ji{x) 



where 7 = | is the ratio of the radii after and before 
shrinking electron disc. The charging energy shows dis- 
tinct forms for different interaction ranges. When the 
interaction range £ is much larger than the radius of 
the disk, £ ^> b, the electron-electron interaction can be 
roughly considered as the bare Coulomb interaction. In 
this case, the electrostatic charging energy is 

E c « 4£- ( 0.290545 - - In |1 - 7 | ) (1 - 7 ) 2 + • ■ ■ (Ry). 

r so V ft J 

In the other limit b 3> £, the interaction is well screened 
and the electrostatic energy has the form 



E c = 4£ 
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where r s and r' s are the inverse density parameters before 
and after the shrinking. The total energy difference can 
be written as 

AE tot = E c {j) + AE g . 

For a small initial radius b <C £, the energy gain AE g 
dominates over the energy loss E Cl thus, there is always 
finite shrinking. 
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FIG. 2. The shrinking distance Ab = b — a versus disc 
radius b for different screening strength, (a) for screened 
Coulomb potential; (b) for Si MOSFETs. The initial den- 
sity parameter r s = 15as. 

However, the above conclusion is not true in general 
as demonstrated in Fig.||. Figure || shows the shrinking 
distance Ab = b—a versus initial radius b. Ab approaches 
a constant for large b. However, the shrinking shown here 
can not be considered as a macroscopic shrinking because 
the typical A6 is only about 2as, which is far smaller 
than the average distance between the electrons, which 
is r s — 15clb for this calculation. Similar behavior has 
also been observed for other values of r s . 

Thus the intrinsic instability is not sufficient to over- 
come the charging energy cost in order to form the elec- 
tron droplet state. The system is in a marginally stable 
situation. However, in real systems, there are always 
disorders. The low frequency component of a disorder 
potential forms the potential landscape, and electrons 
tend to occupy the low potential valleys. We assume 
that around each local minimum, the disorder potential 
is isotropic and slowly varying. We expand the disorder 
potential around the local minimum up to the quadratic 
term. Therefore, we adopt the following simple model 
for disorder potential 

W(r) = V ^, 

where Vq is the potential depth from center to edge of 
the disc. The energy gain by shrinking to a radius a can 
be evaluated as 



AE W = Al J W(r)p(r)dS 

= y (t 2 - 1) ■ 

The total energy difference in the limit b 3> £ can be 
written as 

^- = -(y -4r) d-7 2 ) + A^. 

The effect of the electrostatic energy will be suppressed 
by the potential fluctuation. Large value of Vq is hi favor 
for the disc to shrink. 

AE tot (Ry) Ab (a B ) 
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FIG. 3. (a) The shrinking distance A6 = b — a versus 
radius b for Si MOSFETs with D = 10ag. Top curves are 
for larger disorder potential Vb's. (b) AEtot versus radius b. 
Lower curves are for larger disorder potential Vb's. The initial 
electron density parameter r s = 15as. 

Fig. 3 (a) plots the Ab = b — a as a function of the initial 
radius b for Si MOSFETs with D = 10a B and r s = 15a B - 
There exists a critical Vq (~ 0.08) above which Ab cx b. 
Thus, for large b there is a macroscopic shrinking for 
Vo > V C . Similar result has also been obtained for the 
screened Coulomb interaction with 1/A places the role of 
D. In Fig. 3(b) we plot the energy change as a function 
of the initial radius b. It is clear that larger value of Vq 
gives rise to larger energy gain. We have carried out the 
calculations for many values of r s and the resulting phase 
diagrams are plotted in Figure 4. 

Figure |] shows the phase diagram in the r s — £ plot for 
the case of b ^S> £. Each value of Vo corresponds to a curve 
in the figure. The curves for larger Vq are above those for 
smaller Vo- On the right side of the curve for a given Vo, 
the electron disc will have a macroscopic shrinking, thus, 
an electron droplet phase is stable. To form the electron 
droplet state, the screening of the electron-electron in- 
teraction and the potential fluctuation are both crucial. 
The smaller interaction range between the electrons and 
the lower electron density, the easier to form the droplet 
phase. 
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FIG. 4. The phase diagram in the r s — £ plot. r s is the 
initial density parameter of the electron disc. D or 1/A rep- 
resents the interacting range £ of the electron-electron inter- 
action. Each value of Vo corresponds to a curve in the figure. 
On the right side of the curve, the electron disc will have a 
macroscopic shrinking. 

Before summarizing, we would like to make several 
comments, (i) In this paper, we only consider one elec- 
tron disc, corresponding to one drop in the droplet state. 
In real systems, electrons tend to occupy valleys of poten- 
tial fluctuations to give rise many such drops. The size of 
each drop is determined by local potential depth, (ii) We 
only consider zero temperature effect in this paper such 
that the "gas" phase is empty. At a finite temperature, 
The "gas" phase is occupied by lower density electrons. 
Thus, a finite temperature enhances the possibility for 
the droplet state since the density difference between the 
liquid and the " gas" is smaller and the charging energy is 
less costly. However, in order to form the droplet state, 
the temperature has to be below the cohesive energy (the 
energy cost to remove an electron from the liquid phase 
to the "gas" phase [Q). (hi) We believe that the recently 
observed 2D metal-insulator transition might be the per- 
colation transition of the liquid phase in the droplet state 
Pj . The percolation here is semi-quantum in nature, dif- 
ferent from the conventional classic percolation [|l0| . (iv) 
In order to have a percolation, the "gas" phase needs to 
have a much smaller local conductivity than the liquid 
phase. This requires that a typical length scale of the 
" gas" region is larger than the localization length of the 
"gas" phase. The "gas" phase is low in electron den- 



sity which gives rise to a short localization length. Thus, 
one may not need a large shrinking of electron drops to 
realize a percolation transition. We suspect that in ex- 
perimental samples, a typical drop size is of the order of 
fim. 

In summary, we have demonstrated that it is possible 
to have a droplet phase for 2D electron systems at low 
densities. Both electron-electron interaction and disorder 
potential fluctuations are important for the formation of 
the droplet phase. 
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